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Abstract 

We elucidate a non-conserved relaxational nonequilibrium dynamics of a 0(2) symmetric model. 
We drive the system out of equilibrium by introducing a non-zero noise cross-correlation of ampli- 
tude Dx in a stochastic Langevin description of the system, while maintaining the 0(2) symmetry 
of the order parameter space. By performing dynamic renormalization group calculations in a 
field-theoretic set up, we analyze the ensuing nonequilibrium steady states and evaluate the scal- 
ing exponents near the critical point, which now depend explicitly on ■ Since the latter remains 
unrenormalized, we obtain universality classes varying continuously with ■ More interestingly, 
by changing Dx continuously from zero, we can make our system move away from its equilibrium 
behavior (i.e., when Dx = 0) continuously and incrementally. 
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I. INTRODUCTION 



The concept of universality near critical points in equilibrium systems has a long history 
and is theoretically well-developed [ij]. When equilibrium systems undergo second order 
phase transition at a critical point, they display universal scaling properties for thermo- 
dynamic quantities and correlation functions. These are characterized by a set of scaling 
exponents, which are universal in the sense that they depend only on the spatial dimen- 
sion d and the symmetry of the order parameter (e.g., Ising, XY etc.) pj], but not on the 
parameters that specify the (bare) Hamiltonian. Notable exceptions are the 2d XY model 
and the related models, where the renormalization group flow is characterized by a fixed 
line and consequently the scaling exponents exhibit a continuous dependence on the value 
of the bare stiffness parameter that appears in the model Hamiltonian. The idea of univer- 
sality may be readily extended to equilibrium dynamics close to critical points, where the 
systems exhibit universality through the dynamic scaling exponents, which characterize the 
time-dependent unequal-time correlation functions. Their universality classes depend upon 
the presence or absence of conservation laws and the non-dissipative (reactive) terms in the 
underlying dynamical equations ^j. For driven dissipative out of equilibrium systems with 
nonequilibrium steady states (NESS), the general picture about universality is still wide 
open. In the recent past, attempts with significant success have been made in classifying 
the physics of non-equilibrium systems at long time and large length scales into universality 
classes. For example, the robustness of the standard universality classes in critical dynamics 
to detailed-balance violating perturbations are shown in Ref. 3|. In addition only models 
having conserved order parameters and spatially anisotropic noise correlations exhibit novel 
features. In contrast, recent works demonstrate that truly non-equilibrium dynamic phe- 
nomena, whose steady states cannot be described in terms of Gibbsian distributions, are 
rather sensitive to all kinds of perturbations. Well-known exarnples include driven diffusive 
models 4], and fluid- and magnetohydrodynamic- turbulence 5|-[3]. Overall, in contrast to 
equilibrium systems, how one may classify the universality classes for systems out of equi- 
librium remains an unresolved issue till the date. It is well-known that for models driven 
out of equilibrium, not only dynamical properties but even the static properties (e.g., the 
static correlation functions) depend crucially on the distributions of noises which appear in a 
Langevin description of the model. In light of this a useful strategy is to investigate nonequi- 



librium universality classes is to construct simple models with non-thermal noises, whose 
dynamics will reveal this sensitive dependence of universal properties on noise distributions 
in a systematic manner. 

In this article we examine the particular issue of universality in non-equilibrium for a 
simple relaxational dynamics (model A in the terminology of Ref. {2 1) of an 0(2)-symmetric 
system (equivalently, the classical XY model) in dimension d = dc — e where the upper critical 
dimension of the model dc = 4. The 0(2) symmetric model is a special case with N = 2 
for the more general 0{N) model. The equilibrium critical dynamics of these models are 
discussed in details in Refs. The model is forced out of equilibrium by specific choices 
of the variances of the additive noises in the Langevin equations for the dynamical variables 
(see below). Phase transitions and associated universal properties at the critical point in 
systems with relaxational (model- A type in the language of Ref. |2]) dynamics have been 
shown to be remarkably robust against the introduction of various competing dynamics 
which are local and do not conserve the order parameter Q], including those which breaks 
the discrete symmetry of the system 9|. We show that its NESS depend sensitively on 
parameters of the model. We use field theoretic renormalization group calculation 
using dimensional regularization [l^ based on an e-expansion 11 1 scheme. 
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Our principal results are: (i) As a temperature-like variable in the model (see below) 
is lowered, our model undergoes a phase transition from a high temperature paramagnetic 
disordered phase to a low temperature ferromagnetic ordered phase undergoing a second 
order phase transition at a nonequilibrium critical point, (ii) Universal scaling behavior 
near the critical point determined by a set of standard scaling exponents characterizing the 
correlation and the response functions that depend explicitly on the magnitude of the noise 
cross- correlations; in effect we obtain a continuous universality parameterized by the noise 
cross- correlations, the latter being a marginal operator in the model. The remainder of the 
paper is organized as follows: In Sec. |Tl]we set up our continuum 0(2) symmetric dynamical 
model for a non-conserved order parameter to study its universal properties near the critical 



point. We introduce noises which break the Fluctuation-Dissipation-Theorem (FDT) |15| . 
and thus drive the system out of equilibrium, but keep the rotational invariance in the order 
parameter space unbroken. In the next Sec. IIIII we set up the field theoretic formulation 
for our model in terms of a path integral description. We use a diagrammatic perturbation 
theory and calculate fluctuations corrections to different vertex functions up to the two-loop 



order. We then use a minimal subtraction scheme to calculate different critical exponents 
within an e-expansion. In Sec. |IV]we summarize and discuss the implications of our results. 



II. MODEL EQUATIONS 

In this Section we set up our model equations to describe a simple nonequilibrium gener- 
alization of the relaxational (model A) dynamics in the overdamped limit of a non-conserved 
0(2) symmetric order parameter. The equilibrium characteristics of this dynamical model 
have been extensively discussed in the literature, see, e.g., Ref. j2|. We consider a second 
order phase transition described by a vector order parameter 0j, z = 1, 2. As we further- 
more assume isotropy in order parameter space, the static critical properties are described 
by an 0(2)-symmetric 0^-type Landau- Ginzburg- Wilson free energy functional in d space 
dimensions, 

m] = J d'x[^{<Pr' + 02^) + \{{V<P,f + (V02)n + + (1) 

where r = {T — Tc)/Tc is the bare relative distance from the mean-field critical temperature 
Tc and m > is a (bare) coupling constant. The free energy functional F is manifestly 
rotation invariant in the order parameter space. This F determines the equilibrium prob- 
ability distribution for 0j. Free energy functional F allows us to compute any of the two 
(independent) critical exponents, e.g., the anomalous dimension rj and the correlation length 
exponent u, by means of renormalization group procedure, based on diagrammatic pertur- 
bation theory with respect to the static non-linear coupling u within a systematic expansion 
in terms of e = 4 — about the static upper critical dimension dc = 4. These exponents 
have well-defined physical meaning. For example, the exponent t] characterizes how the 
order parameter correlation function at criticality decays in a spatial power-law fashion, 
(0i(r)0j(r')) oc l/|r — r'\'^~'^~^^6ij, or equivalently of the static susceptibility x(ct) oc \/q^~'^ 
where q is a wavevector, and the exponent v describes how the correlation length ^ diverges 
as the renormalized critical temperature Tc is approached, ^ oc |T — Tc\^^ . Further, the 
fluctuation-corrected true transition temperature is smaller as compared to the mean- 
field critical temperature Tqc , i.e., Tqc = Tc — Tqc < 0. 

In contrast to equilibrium systems, for systems out of equilibrium, there is no detailed 
balance and even the static quantities must be calculated from the underlying dynamics 



directly. The description of the dynamics of such systems in terms of continuum degrees 
of freedom are often based on stochastically driven Langevin equations of motion for the 
relevant dynamical degrees of freedom. For Langevin equations describing processes relaxing 
towards a thermal equilibrium state the correlation functions of a given degree of freedom 
and the corresponding susceptibility are connected through the FDT which in turn fixes 
specific relations between the variances of the noises and the diffusivities. For example, the 
non-conserved relaxational (model A) dynamics for a vector order parameter (pi is given by 

where z = 1,2; F is given by ([T]), F is a kinetic coefficient and Qi are temporally delta- 
correlated zero-mean Gaussian stochastic noises with specified variances. Assuming spatial 
translational invariance we can write generally 

{g^{cl,t)g,{-q,0)) = 2Di,{q)6{t). (3) 

If we now set Djj(k, t) = 2KBTr6ij, where Kb is the Boltzmann constant and T is the 
temperature, then the FDT is obeyed and the corresponding Fokker-Planck equation admits 
a steady-state equilibrium solution P^q ~ exp[—F/ KbT]. In contrast, in nonequilibrium 
situations there are no general relations linking the noise variance and the kinetic coefficients 
and the FDT is broken. Since noises in a Langevin description describe the effects of the 
environment (e.g., thermal baths), such nonequilibrium noises reflect external drives. What 
are the simplest choices of the noise variances which explicitly break the FDT, without 
having to break the 0(2) symmetry? One possible way to do that is to introduce two 
different noise strengths in the noise correlation matrix and break the FDT. This can be 
realized by the choice 

{gi{ci,t)gii-q,0)) = 2rD,6{t), 
(52(q,t)^72(-q,0)) = 2rD25(t), 

(^7i(q,t)^72(-q,0)) = 0. (4) 

Such a choice as above will certainly break the FDT but unfortunately will break the 0(2) 
symmetry of the ensuing dynamics as well. Ref. |l6| investigated nonequilibrium critical 
properties of 0(n)-symmetric models with reversible mode-coupling terms. Speciflcally, a 
variant of the model of Sasvari, Schwabl, and Szepfalusy (SSS) is studied, where violation 



of detailed balance is incorporated by allowing the order parameter and the dynamically 
coupled conserved quantities to be governed by heat baths of different temperatures. They 
however find that upon approaching the critical point detailed balance is restored, and 
the equilibrium static and dynamic critical properties are recovered. Yet another option is 
to couple the system with the corresponding conserved angular momentum and introduce 
dynamical anisotropy in the noise for the conserved quantities, i.e., by constraining their 
diffusive dynamics to be at different temperatures T|j and T± in (iy- and c/^-dimensional 
subspaces, respectively, see Ref. for detailed calculation for the SSS model for planar 
ferro- and isotropic antiferromagnets. Ref. showed that the equilibrium fixed point 
(with isotropic noise) to be stable with respect to these non-equilibrium perturbations, and 
the familiar equilibrium exponents therefore describe the asymptotic static and dynamic 
critical behavior. Novel critical features are only found in extreme limits, where the ratio 
of the effective noise temperatures T\\/Tj_ is either zero or infinite. In a similar study, 
Ref. [l8| discussed nonequlibrium dynamics in a liquid-gas model with reversible mode 
couplings. The model is driven out of equilibrium by introducing different temperatures 
for different dynamical variables, or, by having anisotropic noises. However, no new genuine 
nonequilibrium stable fixed point is found (within one- loop calculations). Similar approaches 
to nonequilibrium critical dynamics of the relaxational models C and D (in the language of 



Ref. [2j]) are discussed in Ref. 19| and involve coupling a non-conserved and conserved order 
parameter, respectively, to a conserved density, with the order parameter and density fields 
are being in contact with heat baths at different temperatures. Within a one-loop calculation 
it finds, in certain cases, continuously varying static and dynamic critical exponents, as a 
function of a dimensionless nonequilibrium parameter in the model. An alternative route 
to introduce detailed balance violation in the simple model A-type relaxational dynamics 
for the order parameter 0jis to introduce non-zero noise cross correlations which will make 
the noise matrix off-diagonal. This will break the FDT as the noise matrix is then not 
proportional to the kinetic coefficient matrix (which is proportional to the unit matrix in 
the present case). We take cross noise strengths as -D(q). We write 

(^i(q,i)^7i(-q,0)) = (^2(q,t)(72(-q,0)) = 2OT5(t) 

{gi{q,t)g2{-q,0)) = 2D{q)r6{t). (5) 
In general the function D[q) is a complex function of wavevector q. 



The form of the function -D(q) may be restricted by demanding rotational invariance 
of the noise variance matrix (equivalently by demanding 0(2) symmetry of the dynamics). 
Under a rotation by an arbitrary angle 9 in the order parameter space the noise variance 
matrix transforms to 

cos 6 sin 6 \ I D D \ I cos — sin , 

^' = r . „ 1 U I . „ „ ■ (6) 




— sin6' cos6' / \ D* D j \ smO cos 9 
where the noise variance matrix before rotation is 




(7) 



Now we demand N = N' due to rotational invariance. This after a simple algebra then yields 
that the noise cross correlation amplitude should be fully imaginary or -D(q) = —D*{q). 
Since in the real space -D(r) must be a real function, we find -D(q) must be an odd function 
of q. In order for the noise cross correlation to have the same naive dimension as D (so that 
both D and D are equally relevant in an RG sense), we set D{q)D{q) = D\ where D\ is a 
constant (and has the same dimension as D'^). We henceforth replace D{(\) by i-D(q) where 
D(q) is now completely real. This reflects the fully imaginary nature of the cross correlation 
explicitly. Thus the explicit forms of the two equations of motion for 0i and 02 are 

f~dr ^ '^'^^ + CV 01 - -0f - -0102 + 

1 502 , ^9 , W , o M , , 9 92 / N 

f ^ = -^^2 + CVV2 - ^0^ - ^020? + ^, (8) 

complemented by the noise variances as below: 

(^i(q,t)^?i(-q,0)) = (^?2(q,t)^?2(-q,0)) = 2DT5{t) 

(^7i(q,t)^?2(-q,0)) = 2z^(q)r5(t). (9) 

One may in addition consider including a conserved angular momentum as a slow variable 
in the problem (see, e.g., model E in Ref. {2!). We do not do that here for simplicity. 
Model equations ([H]) suffices for our purposes of exploring nonuniversal features in a simple 
setting. Are there limits on the value of in this model? To obtain that we demand the 
noise variance matrix to have eigenvalues which are real positive or zero. The eigenvalues 
concerned are D ± Dx- Thus \Dx\ < D, or in terms of a dimensionless number A^x = 



(Dx/D)"^, Nx < 1. In the subsequent calculations we will find that A^x enters into the 
expressions of different scaling exponents explicitly. 

Equations of motion are written in an 0(2) invariant representation. Using the 
equivalence between 0(2) and U{1) representations, one may write an equivalent U{1) rep- 
resentation of the dynamics. The free energy in the U{1) representation takes the form 



Fuii^r] = J d^xlrtpr + (VV^)(V^/^*) + ^{^r)% (10) 

where complex fields ip = -^{(f)i+i(f)2)] V'* is the complex conjugate of tp. The corresponding 
Langevin equations of motion in the overdamped limit are given by 

where zero-mean Gaussian distributed complex noise ^ has the following correlations in the 
Fourier space: 

(e(q,t)^(-q,0)) = = (r(q,t)r(-q,0), 
(e(q, i)r(-q, 0)) = 2DT5it) + 2z^(q)r5(t). (12) 

Thus introduction of noise cross-correlations in the 0(2) description is equivalent to adding 
an imaginary and odd function of q in the variance (^^*). Before we embark upon detailed 
calculation let us consider possible physical (microscopic) realizations of our continuum 
model in terms of stochastic lattice-gas models. However, what we discuss below does not 
fully and precisely define a microscopic model, but rather outlines broad features that an 
eventual appropriate microscopic realization should posses. Consider a system of XY ((0(2) 
spins) either on a (hypercubic) lattice or a continuum in c/- dimensions, interacting with 
an additional mobile species in the system which diffuses randomly, undergoing symmetric 
exclusion process (SEP) to any of the nearest sites, if vacant. A simple model of interaction 
between these two species could be where each diffusing particle carries an XY spin attached 
to it, and the nearest-neighbor exchange coupling Jij that defines the XY model is related 
to the local particle density ni{t) at site i via Jij oc ni{t)nj{t). Next, noting that the 
microscopic dynamics of both the spins and particles are stochastic, characterized by two 
sets of random numbers gii{t) and g2iit), respectively, we impose that gi and g2 are cross- 
correlated, with the cross-correlation function being of the form (in the continuum limit) 
y45(xi — X2) + -B(xi — X2), where xi and X2 are two points in the lattice, A is a numerical 
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constant and i?(x) is an odd function of position x having the same dimension as 5(xi — 
X2). The presence of the odd function -B(xi — X2) in the cross-correlation function ensures 
lack of reflection invariance of the underlying stochastic microscopic dynamics. Thus the 
measured quantities (e.g., correlation functions of appropriate densities) should reflect this 
lack of reflection invariance. At this level, the dynamics for the additional species is clearly 
conserving. This implies there will be no timescale present on which the additional variables 
of the mobile species can be treated as fast and eliminated to yield an effective equation of 
motion for the XY spins alone with the effects of the diffusing species buried in the additive 
noises in the effective spin equations. However, if particle nonconservation is introduced, 
e.g., via evaporation-deposition affects, the local particle density dynamics will be fast and 
then may be eliminated to produce an effective spin dynamics. Since the noises in such 
effective theories contain information about the already eliminated fast degrees of freedom (in 
this case the local diffusing particle density), there will be non-zero noise cross-correlations 
of specified structures as above, due to the particular chosen structure of the underlying 
reflection invariance breaking microscopic dynamics. Alternatively, one may introduce the 
driving as a temporally delta-correlated fluctuating magnetic field h(x, t) = {hx,hy) with 
hx and hy having short ranged spatial correlations as in iQ. In both the cases, noise cross- 
correlations of appropriate structures will be generated in the effective Langevin description. 
With this short background in mind, let us now investigate the universal scaling properties 
of the model described by the Langevin equations (JHl) together with the noise variances (l9l). 
The presence of non-linear terms in Eqs. (IHl) rules out exact solutions, and we resort to 
perturbative calculations that we discuss below. 

III. NONEQUILIBRIUM STEADY STATES 

Let us first consider the high temperature phase of the system. At high temperature with 
r > 0, the system is in the paramagnetic phase, i.e., (0j(x, t) = 0, where (...) means averaging 
over the noise distributions. The correlation length ^ remains finite for all r > 0. The only 
effect of the noise cross-correlations is to make the cross-correlation function (0i(x, t)(f)2{0, 0)) 
non-zero with a finite correlation length ^. Further, as in equilibrium critical dynamics, 
the paramagnetic phase is linearly stable and the fluctuations have a finite life time for 
all wavevector. Nevertheless, the FDT is violated for all r > due to the noise cross- 
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correlations. 

Near the critical point, the system becomes scale invariant and the correlation length ^ 
diverges, leading to the emerging macroscopic physics near the critical point being vastly dif- 
ferent from the paramagnetic phase. A quantitative description of the nature of correlations 
near the critical point requires the principles and formalisms of the Dynamic Renormal- 
ization Group (DRG), which we execute here by using a field-theoretic framework. Detail 



discussions of the technical aspect o: 
in the literature, see, e.g., Refs. [10 



field-theoretic DRG calculations are well-documented 



20|. In order to set up the background let us examine 
the linearized version of the model equations ([H]) together with the noise correlations (EJ at 
r = by dropping all the non-linear terms {u = 0). The system remains 0(2) invariant, 
but the FDT is already broken at this linear level due to the noise cross-correlations. Obvi- 
ously, the field correlations from the linearized model equations can be exactly calculated. 
In this linear theory, in the critical region, defined by r = 0, the linear theory is massless 
resulting in divergent long wavelength fluctuations, as can be seen by explicit calculations 
of the correlation functions Cij = (0j(x, )f:)0j(O, 0)), i = 1,2, which may be written down in 
a scaling form at the critical point r = 0: 

Cn(x,t) = C22(x,t) = x'-^'fit/xn, (13) 

where / is an analytic function of its argument. For the cross-correlation function Cu (and 
by symmetry C21) displays the same scaling form, but with a different amplitude and is an 
odd function of x. 

What is the nature of these diverging fluctuations when the non-linear terms are present 
(m > 0)? The presence of nonlinear terms no longer allows for exact solutions, in contrast to 
the linearized theory. However, this question may be systematically addressed via standard 
implementation of DRG procedure, based on a perturbative expansion in the small coupling 
u about the linear theory. The perturbative corrections to the correlation function may 
be equivalently viewed as arising from modifications (renormalization) of the parameters 
r, u, r, D and the dynamical fields (pi and 02- Renormalizability of the theory ensures that 
correlator Cij will retain scaling forms similar to f|T3l) with exponents different from those 
appearing in (fT3l) and new scaling functions at (renormalized) r = 0: 

Cn(x,t) = C22(x,t) = x^-''~^Mt/x^), (14) 



10 



where t] and z are the anomalous dimension and dynamic exponents respectively, and fs 



is a new scaling function 2l|. For the linear theory described above, 1] = and 
The nonlinear coupling u is expected to change these exponents for the linear theory. The 
equilibrium critical dynamics of several nonlinear problems have been described in Ref. j2|. 
In our subsequent analysis below we assume renormalizability and justify it post facto by a 
low order (up to two-loop) perturbation theory. 

Operationally, the DRG procedure is conveniently performed in terms of a path integral 
description based on a dynamic generating functional which is to be constructed out of the 
Langevin equations (Ml and the corresponding noise variances given by Eq. following 
standard procedures |l3, 22|. The dynamic generating functional for the present model is 
given by 

f ^1 ^1 r D [ d^k D f d'^k ^ . r d'^k , D 

Z= / T>(piV(p2T>(piT>(p2exp[-— / 77r-T70i0i " / 77^-^70202 



r J {27rY'^'^' r J {2ttY^"^' "J (27r)^^' r 
f d'^k 2 A ^ ^ 5F^^ f d^k 2 rl ^ SF., 

^ J ^^'^^ + W.^^ J J^'^^'^f^'"^' + 50^^]' ^^^^ 

where 0i and 02 are conjugate response fields corresponding to the order parameter fields 0i 
and 02 used to obtain the noise-averaged generating functional Z. Clearly, with our choice 
of the noise cross-correlations, the generating functional respects the 0(2) symmetry of the 
underlying dynamics, but explicitly breaks the FDT. It now remains to be seen whether 
this breakdown of the FDT remains valid even for the effective, renormalized version of this 
theory, or it ^restored in the long wavelength limit. In the above we have used the Ito 
prescription [23] while writing down the generating functional ( IT5|) . 

The perturbative calculational framework begins with the construction of the perturba- 
tion expansion of one-particle irreducible Feynman diagrams for all possible correlation and 
the corresponding vertex functions constructed out of the fields 0i, 0i, 02 and 02. Such 
a perturbation expansion is meaningful only wheen the coupling u is small. This is gen- 
erally accomplished using the standarad DRG procedure which involves an order by order 
expansion about d^. — rf, where d^ is the upper critical dimension and d the physical space di- 
mension. As a result, the renormalized coupling ur flows to small values, of order e = dc — d. 
As for equilibrium field theories, a straight forward scaling analysis yields that dc = 4: for 
this model: We scale x — > b^x and t — )■ btt in the action, where bx > 1 and bt > 1 are 
arbitrary parameters, and find out how the other quantities scale to maintain scale invari- 
ance. It is seen from the bare equation leaving aside the non-linear terms we must have 
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bt — bl — 6^ in order to have dynamical scaling, where we have taken bx — b. For the action 
to remain invariant, the fields (pi and 02 pick up a canonical dimension d/2 — 1 and the 
coupling constant u a dimension A — d. Hence the critical dimension at which the coupling 
constant u becomes dimensionless is dc = 4. This yields that the perturbation theory will 
be infra-red (IR)-singular for d < 4, and consequently the system will show non-trivial crit- 
ical behavior in that regime, while for d > 4 the perturbation theory contains ultra-violet 
(UV) divergences, and the (static) mean-field exponents together with dynamic exponent 
z — 2 (dynamic exponent of the linearized theory) will describe the system at the critical 
point. To make the field theory UV renormalized it is needed to introduce the multiplica- 
tive rcnormalization constants in order to render all the non-vanishing two- and four-point 
functions finite. Within the DRG procedure this is achieved by demanding the renormalized 
vertex functions in the theory, or their appropriate momentum and frequency derivatives, 
to be finite when the corresponding loop-integrals representing fiuctuation corrections are 
considered as functions of conveniently chosen frequency and momentum, well outside the 
IR regime. In order to compute the associated one- and two-loop momentum integrals we 
employ the dimensional regularization scheme and choose r = /x^ as our normalization point, 
where /j, is an intrinsic momentum scale of the renormalized theory. From the rcnormaliza- 
tion constants (Z-factors) that render the underlying field theory finite in the ultraviolet 
(UV), one may then derive the RG fiow equation, which describes how correlation functions 
change under scale transformations. Since the theory becomes scale-invariant in the vicinity 
of a critical point (or an RG fixed point), one may employ the information previously gained 
about the UV behavior to access the physically interesting power laws governing the infrared 
(IR) regime at the critical point ( r oc T — — )■ 0) for long wavelengths (wavevector q — )■ 0) 
and low frequencies {u — )■ 0). The scaling behavior of the correlation or vertex functions 
may be extracted by finding their dependence on /i by using the RG equation. 

Ward identities due to the rotational invariance of the model (in the order parameter 
space) ensures the following exact relations between different vertex functions: 



r0i,^i(k,a;) = r^^^^{k, uj) = rii(k,a;), 

r^.0,(k,^) = r^,^,(k,u;) = r2o(k,c^) 



(16) 
(17) 
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and 

r0i0i^i^i(ki,k2,k3, -ki -ka -k3,a;i,a;2,c^3,-c^i -c^2 -ws) 

^ r02<^2<^2<^2('^l' ^2, kg, -ki - ka - kg, UJi,UJ2, U3, -Ui -UJ2- Us) 
= 2r^i,^^^^^^(ki,k2,k3,-ki -k2-'k3,UJi,UJ2,UJ3,-UJi-UJ2-UJ3) 

^ 2r^^^^^^^^(ki,k2,k3, -ki - ka - kg, wi, ^2, ws, -wi -^2-6^3) 

= ri3(ki, ka, ks, -ki - ka - kg, wi, a;2, o-s, -uji -U2- u^) (18) 

Thus in the present model, the only UV-divergent two- and four-point vertex func- 
tions which require multiplicative renormalization are (i) (9(^rii(k, cu), (ii) 5fc2rii(k, cu), (iii) 
rii(k,cj), (iv) r2o(k,cj) and (v) ri3(ki, ka, k3, -ki - ka - k3, wi, ^2, tua, -^^i -^2-^3). 
Each of them is to be rendered finite through multiplicative renormalization by means of 
introducing a Z-factor. Thus there are 5 Z-factors in total. However, there are four pa- 
rameters (r, D, T, u and two fields (0^, 0^, i — lor, 2), thus six altogether, available for 
renormalization; thus this leaves us at liberty to choose one of the renormalization constants 
in a convenient manner. 

The renormalized kinetic coefficient F/j, noise strength D^, mass Tr and coupling constant 
ur are defined in terms of the above vertex functions as 

9<.rn(o,o) = -^, 

J- R 

2Dr 

dk^TiiiO, 0) = 1, rn(0, 0) = TR, r2o(0, 0) = — ri3(ki = O, a;, = O) = UR. (19) 

i R 

The above definitions of the renormalized parameters allow us to calculate different renor- 
malization Z-factors in the problem. 

The perturbation theory here is constructed out of the bare propagator and correlation 
functions which are to be read off from the harmonic part of the action functional. From 
the generating functional we get the bare propagators as 

(0i(k, co)M-K -co))o = G'lik, u) = ^ ' = (02(k, co)M-K -^)) = Gl{k, u) 

P — r T -r fv 

(20) 

and bare correlators as 

2D 

i-[f2 + [T + ) I 

{Mk,u^)M-k, -a;))o = C°(^,a;) = ^^^^^^ . (21) 

i-[Y2 + [T + K ) \ 
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The self energy T,a(k.,u!) is formally given by the Dyson equation: 



G^\k,uj) = -'-^ + T + e - J:G{k,co) = G^\k,uj) = Ti^{k,uj). (22) 

In the same way one may define I^]:){k,uj) through the relation 

r2o{k,uj) =2D + J:D{k,uj). (23) 

We now calculate fiuctuation corrections to the relevant vertex functions. One-loop dia- 
grammatic contributions to Fis do not receive any contribution from Dy^ and are structurally 
identical to their equilibrium counterparts. Similarly one- loop corrections to E(3(k, a;) are 
independent of Dx and are identical to the corresponding equilibrium contributions. In 
contrast, there are additional two-loop Dx -dependent diagrammatic corrections to Eg (A;, a;) 
and Sd(/c,C(;) whose evaluations require careful consideration. 




FIG. 1. Representative two-loop diagrams coming from non zero noise cross correlations contribut- 
ing to 'EG{k,oj) (left) and 'ED{k,oo) (right). A line with a filled circle represents a correlation 
function, a line without any filled circle represents a propagator. We do not show all the diagrams 
here. 

We separately consider the following contributions di^'EaiO, 0), dk2'EQ{0, 0) and Ed(0, 0). 
The cross-correlations contributions 'EQ{k,u!) and 'E^{k,uj) to 'EaikjOu) and 'EoikjUj) re- 
spectively, which do not arise in equilibrium, are all even in D{q} (or in D^), since the 
model must be invariant under •<-)■ — -Dx- Such contributions to EG'(k, cj) are of the form 
(up to numerical factor): 

d% d% Z)(qi) Z)(q2) r 



=sft")=«^(^i)/ 



{27tY {27tY (r + qf) (r + ql) [-tu + Tqj + Tql + T{3t + (k - qi - qa)^}] 

(24) 
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and 



{27iY {27iY (r + qf) (r + q^) r + (qi + qs)^ 

(25) 



3r + g? + gi + (qi + q2)2' 
where k and u are external wavevector and frequencies, respectively. We separately 
need to find out the k^u and k'^u^ parts of of the integral We need to calculate 

-^T,^(k, u)\k=o,ui=o where a = k'^, u. Let us consider a = u: 

d f d% d% D{q,) D(q2) 1 

(26) 



dco J {27rY {27iY (r + qf) (r + q^) -ico + T[3t + qf + g| + (qi + qs)^] ' 
Since we are using an e-expansion based on a minimal subtraction scheme, we need to extract 
the diverging parts of ( l26l) . to be given by poles in e. It is noteworthy that integral ( l26ll 
has a structure very similar to and has the same logarithmic divergence (by simple power 
counting) as its equilibrium counterpart [i.e., when -D(q) is replaced by D in ( |26l) ]. Clearly, 
the dominant contribution to it comes from p = qi + q2 ~ 0, which controls the critical 
behavior of this integral. We write (retaining only the small-p contribution), up to constants 
and numerical factors 
5^x„ M 9 f d% d% D(qi) D(p-qi) 1 



^^Eg(A;,u;)|fc_o,^_o ^ J ^2nY {2ttY (r + qj) (r + qi) -luo + r[3r + ql + ql + f] 

d_ r d% d% Djqi) D(qi) 1 

^ dcu J i2nY {2nY (r + qj) (r + q^) -icu + T[3t + qf + q^ + 

= n2 A /" ^'^^i ^"^^2 1 1 1 



X 



du J {2tiY {2'kY (r + ql) (r + ql) -icu + T[3t + ql + q^ + (qi + qa)^] 
Dl f d'^qi d'^q2 111 1 



-I ^ 



312 J ^2tiY (27r)'^ (r + g?) (r + gf) r + (qi + q^)^ 3r + + ql + (q, + q2)2 ' 

(27) 

where, we have replaced integral fl2B]) by its dominant contribution coming from p ~ 0. The 
last line of Eq.f l27p is obtained by symmetrizing the previous line. The reduced integral 
( 127|) may now be evaluated exactly in the same way just as its equilibrium counterpart (i.e., 
when is replaced by D). In a similar way, cross-correlation contribution to the anomalous 
dimension (again logarithmically divergent on a naive power counting basis) may be written 
as (up to constants and numerical factors) 

•9 M r^2 d f d% d'^q2 1 1 1 



Q}^2 OK ,-7,«-u,.-u xQf^2j (^27fY{27TY{r + ql){r + ql)T + {k + qi + q2)^'' 
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and the cross-correlation contributions to noise strengths (l25l) become (up to constants and 
numerical factors) 

EUO,0) ^ DD^ f '''^ ''^^ 111 



" J {27rY {2nY (r + gf ) (r + qi) r + (qi + qa)^ 3r + qf + ql + (qi + qs)^ ' 

(29) 

In the above, in our evaluations of the cross-correlation contributions ^S^(A;, (X')fc=o,ij=0; ci = 
fc^, oj and S^(0, 0), we have picked up the dominant contribution given by p = qi + q2 ~ 0. 
Subdominant contributions are neglected and are expected to be small as we heuristically 
justify: For example, for p ^ qi, the integrand in S^(0,0) (see Eq. ([29] above) is 

D(qi) ^(q^) D 1 D(qi) D(p) D 1 

(r + ql) (r + ql) t + (qi + qs)^ 3r + ql + ql + (qi + q2)2 (r + g^) r + p2 ^ + p2 37- + 2p2 ' 

(30) 

which can be both +ve and -ve since -D(qi) and -D(p) are odd functions of their arguments, 
and hence contributions from outside the dominant region p ~ will be small due to mutual 
cancelations. Our results, although backed up by heuristic arguments, nevertheless bring 
out remarkable new features, as we shall see below. Thus, after putting every diagrammatic 
contribution (up to two- loop order) together we obtain for SG(k, w) as 

T^cik, u) = — — / ——r — ■ — - + Tu D ' 



X 



3 J (27r)2r + g2 18 9j J {2tt)'^ {2tt)'^ t + ql 

1 1 
T + ql -lUJ + r{3r + ql + g| + (k - qi - qa)^} 

2 1\ /• d^qx d% 1 1 1 



6 9j J (27r)^ (27r)^ r + g^ r + g| -zo; + r{3r + g^ + ql + (k - qi - q^)^} ' 

(31) 



Similarly the two loop contributions to T,d{0,0) comes out to be 



8j J i2n)^ {27r)^ r + gf r + g| 3r + ql + ql + (qi + qa) 



r + (qi + q2)2 T \3 18 J J (271)'^ {2n)'i r + ql r + ql 

X - - (32) 

r + (qi + q2)2 3r + g2 + g2 + (qi + q2)2 ■ 

Although formally there exists two-loop diagrammatic corrections (there are no one-loop 
corrections) to -D(k), all of these vanish in the long wavelength limit due to the fact that 
-D(k) is an odd function of wavevector k. In Fig. ([2]) we consider one such two- loop diagram: 
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FIG. 2. A representative two-loop diagram contributing to Ex(A;,a;). Symbols have meanings as 
before. 

The corresponding expression Il^{k,u) is (up to constants and numerical factors) 

r d\ d% D{qi) 11 1 
'"^^^y {2nY {271Y r + qlT + qlT+ (q^ + q^)^ 3r + ql + g| + (qi + q^)^ " ^^^^ 

Clearly, the integral in (l33ll vanishes. We shall come back to this issue of non-renormalization 
of -Dx(k) again at the end. 

Finally, to complete evaluating diagrammatic corrections we now evaluate the F13 up to 
one-loop order at zero external wavevector and frequency. There are no contributions from 
Z^x to the four point vertex function. We obtain 



F(3 n = u^D ^ ,X {e/2) 



1 1 
4 36 



(34) 



After evaluating all the two point and four point vertex functions we can now renor- 
malize vertex functions F(ii)(0,0), ^F(ii)(ci; = 0,0), g|^F(ii)(0, A; = 0) and F(3^i)(0,0). We 
now define renormalization Z-factors for the present model. We begin by introducing the 
renormalized fields (hf- and chf". 



This implies that the renormalized vertex functions become 



■pi? _ 7-Af/2 7-Af/2p 

^ {N,N) ~ ^ ^ ^ {N,N)- 



We further define the renormalized parameters as 



(35) 



(36) 



= ZdD , = Z^T^i-'' , = ZuuAd^i"-", , F^ = Z^T 



-2 „,R 



d-4 r/J 



(37) 



where /x is the scale factor introduced to make the renormalized parameters dimensionless. 
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Here — — ^— 3- . Thus we get 



d 



1 + . + 



r«o)(o,o) = -i-^ 



9(47r)'^e 27(47r)<^e 



2DV^ 

^ 3(47r)'^e ""^^ ^ 6(47r)'^e 



ln(; 



rf3,i)(o,o) = z-V2z-3/^^ 
rf,i)(o,o) = z-i/2z-^/V 



10^ 

1 - —Duji 
6e 



(38) 
(39) 

(40) 
(41) 



1 - 



3(47r)<^/2e 3(47r)'^e \e 



9(47r)'='e 



+ 1 



= T 



R 



(42) 



from which we can calculate all the Z-factors. We use the freedom to choose one of the 
Z- factors freely to set Zd — 1- Henceforth we set D — 1 for simphcity without any loss of 
generality. Assuming = N^D'^, where A'"x is any dimensionless parameter, the other Z 
factors are obtained up to two loop order as follows 



36 e V 3 / 108 e V 3 

36 e 72 e 3 108 e 

Z = l + —^ ' - — iVv M n - + — A^v 



36 



72 



2(MAd/i^) 1 /7 6\ {uAd^l-'f 



Zr^l 

3 e 18 V2 e 

Defining the Wilson's flow functions as 



108 



54 



7 6^ {uAdii 
2 e 



-t\2 



Cr = /x^lnZr, 
Cd = ii-^lnZo, 

Cr = In — = -2 + In Z^ 



(43) 

(44) 
(45) 
(46) 



(47) 
(48) 
(49) 
(50) 
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and the ^ function for the non-hnear couphng as 



d { 10 



d\x \ 6 

we get a stable nontrivial fixed point at u* = j^e and we can evaluate the critical exponents 
from these flow functions at the fixed point. The flow functions at the fixed point pick up 
values up to order as follows: 









(52) 


— 


50 150 


100 3 


(53) 


Cr = 


-e^ f61n--lV 
50 V 3 / 


—N^e^ f61n--l^ 
150 V 3 ) 


(54) 




-2 + ^e + 0(e2). 
5 




(55) 



The basic renormalization group (RG) equation is derived on the basis that the unrenor- 
malized correlation and vertex functions do not depend on the arbitrary scale fi. Hence if 
we hold the bare parameters D, r and fixed, we must have 

= /x|^[^^/^^^/^r^_^(/.,r^r^x.^)]. (56) 

As Z-factors also depend on the RG equation finally takes the form 



ZN/2zN/2^ 



d ^y, d ^ ^ N d ^ ^ dV^ d dr^ d du^ d 
— + AT/ 2 — In Z + - — In Z + — - — - + -- + 



N,N " 



(57) 



At the critical point we have scale invariance separately under scaling of space, time 
fields and parameters. These are determined by the momentum and frequency canonical 
dimensions of the fields and parameters. After proper scaling as described above in this 
Section we have canonical dimensions of fields and parameters as 

4.0. = - 1 , d;^^^^ = 0, 

d'l , ^ d/2-1 , d'^. . =1 

91,92 ' <Pl,<P2 

4 = -2 , d^^^l 
4 = 0, < = 0. (58) 
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Canonical scale invariance at the fixed point {f3* = 0) for the correlation function implies 

C^{x,t,fi,r^,T^) = and 



a/i ox 



d d 

r 7^ -t—-<K. 

dV^ dt ^ 



R 







(59) 



where ci^ and d!^ are the momentum and frequency dimension of the correlation function 
C(x, t) respectively. In this case = d — 2 and d^ = 0. Therefore from Eq. ( I59l) we have 



r^^c^ibx, b\ r«) = t^c^ix, t, r«) 

The RG equation for the correlation function at the fixed point can be written as 



(60) 



d 



^ ^R d . R d 



C^(x,t,/i,r",r 



R R\ 



(61) 



Combining the two separate spacial and temporal scale invariant equations Eqs. (!59|) and 
using equations f l60|) and fl6T|) we get at the fixed point 



d . R d 

-X^ - CrT^ 



(2 + Cr)t^ + 2-rf-r7 



0, 



(62) 



dx ^ ' Ot^ 

where we have used r] = —Cp. From Eq. (!62|) it can be seen that at the critical point (r = 0) 
and equal time (t = 0) the correlation function should take the form 



Cix) 



~ X 



(63) 



which gives the spatial scaling of the equal-time correlation function at the critical point. In 
case of time- dependent correlation function at the critical point the scale invariant equation 
takes the form 



C^ixAS^) 



0. 



(64) 



Assuming dynamical scaling, the solution of C{x,t) should be of the form C{x,t) 



x^ "^gix?^^^ It\ which implies that 



2 + Cr 



(65) 



should be the dynamic exponent. At equal time (t = 0), near the critical point (r^ ^ 0), 
the equation for C (x, t) can be written as 



c^ixAy 



0. 



(66) 
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This implies that the correlation function should be of the form C{x,t) ~ ''/(x^^/r), 
and the correlation length exponent 



1 



(67) 



For the propagator G 



, the scale invariant equations are given by 
G^{x,t,fi,r^,T^) = and 
G^{x,t,fi,r^,T^) = 0. 



d d 



- 1- 



d' 



G 



dT^ dt 

where (Iq = d — 2 and dQ = 1. From the second of equations fl68|) it is obvious that 



(68) 



d 



l + t 



d_ 

di 



R R\ 



(69) 



The RG equation for the propagator at the fixed point can be written as 



d 1 1 



d 



d 



G^'ix^t^n, r^,r^) 



(70) 



Using equations ( ITOl) and ( ]69l) in ( ]68l) we get 



G'^(x,t,r^,r^) =0, 



(71) 



where 17 = — C^- 

From Eq. f lTTj) the time dependent propagator at the critical point (r = 0) can be written 



as 



-,_|-(2 + Cr){l + 4} + 2-d4„-l,5 



G^{x,t,T^) =0. 



(72) 



Assuming dynamical scaling G{x, t) should be of the form G ~ 3;2-d-q/2-?7/2-2-Cr j rr^+'^r ^ 
Since the static susceptibility x(x) is proportional to Jq°° dtG{x, t) which on integrating over 
time gives us 



(73) 
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From equations (165|) . ( I63i) and (167|) and using equations (IMI) . (152|) and (155|) we get the 
dynamic exponent, anamolous dimension and correlation length exponent of the model to 
the leading order in e: 

ri = —e' + —N,e' + —N^e' In -, (76) 
' 50 150 100 3' ^ ^ 

l = 2-le + 0{e'). (77) 
u 5 

If A^x = we get back the equilibrium exponents as expected. Their contribution from the 
nonequilibrium part depends on the value of A^x- Let us now consider the consequences 
of a non-zero A^x- First of all, as is evident from the results presented above, the static 
susceptibility and the equal-time correlation function do not display the same spatial scaling, 
since t] ^ fj. This is an important evidence of breakdown of the FDT in the renormalized 
theory. Note, in the linearized theory breakdown of the FDT is manifested in the existence 
of non-zero cross-correlation functions: The correlation matrix pick up non-zero off-diagonal 
elements, where as the dynamic susceptibility matrix remains diagonal. Nevertheless, all 
elements of the correlation and the dynamic susceptibility matrix exhibit the same scaling 
properties at the critical point. In contrast, in the renormalized non-linear theory, not only 
the correlation matrix is off-diagonal where as the susceptibility matrix remains diagonal, the 
elements of the correlation matrix scale differently from those of the susceptibility matrix. 
The latter result is purely a non-linear effect. Further, the dynamic exponent for finite 
A^x, z{Nx) is larger than -2(A^x = 0), its value for equilibrium dynamics. Thus relaxation 
for finite A^x is slower than for the corresponding equilibrium dynamics. The correlation 
length exponent u has been calculated only up to 0(e) and is equal to its equilibrium value. 
However, as there are A^x -dependent corrections to Sg(0, 0) at the two- loop order, the value 
of u is likely to be different from its equilibrium value at the two-loop order. Lastly, the 
A^x -dependence of all the scaling exponents are continuous - all of them vary continuously 
with A^x and go over to their equilibrium values when A^x is set to zero. Our claim of 
the scaling exponents varying continuously with A^x rests on the marginality of A^x- We 
have shown explicitly up to two-loop order that -D(k) does not renormalize. Hence A^x does 
not renormalize up to two-loop order and depends on the bare value of . Any non-zero 
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fluctuation corrections to -D(q) must be an odd function of the its wavevector argument. 
In order to have that one must have odd number of internal cross-correlation line. Since 
all internal wavevectors are integrated over, such a contribution will vanish in the limit of 
vanishing external wavevector and frequency. Thus -D(q) and remains unrenormalized and 
hence iVx appears as a dimensionless marginal operator to any order in perturbation. 



IV. SUMMARY 

To summarize, we have analyzed the universal scaling properties of a nonequilibrium 
version of 0(2)-symmetric dynamical model near the critical point. We write down a non- 
conserved relaxational dynamics for the order parameter field. We have introduced cross- 
correlations between the two additive noises in the Langevin equations, so that the FDT 
is immediately broken. We then show that if the cross-correlation is imaginary and odd 
in wavevector, the underlying 0(2) symmetry is still maintained. We calculate the scaling 
exponents of the model in a DRG framework using an e-expansion scheme, where e = 4 — d 
with 4 being the upper critical dimension of the model. We show that at the two-loop 
order there are diagrammatic corrections to the various two-point vertex functions in the 
model arising from the cross-correlations. We have used heuristic arguments to extract the 
dominant contributions to the two-loop diagrams involving cross-correlations, which have 
allowed us to evaluate the respective cross-correlation contributions in a simple and con- 
trolled manner. We finally argue that the cross-correlation amplitude appears as a marginal 
operator in the problem. Since this amplitude appears in the expressions of the scaling expo- 
nents we have an example of a continuously varying universality class. Technically speaking 
we obtain a fixed line, parametrized by the value of the parameter A^x introduced above, 
instead of a single or isolated fixed points. Every point on the fixed line characterizes a 
universality class, parametrized again by N^.. The fixed line begins from A^x = 0, which is 
the equilibrium fixed point. This stands in contrast to, e.g., Ref. jsl, where nonequilibrium 
noises lead to additional fixed points, but not a fixed line as here. There are other dynamical 
models where cross-correlated noises lead to universal properties varying continuously with 
the amplitude of the noise cross-correlations. Notable examples are the stochastically driven 



generalized Burgers model (GBM) j2J| and magnetohydrodynamic turbulence (MHD) |25 |. 



However, the ci-dimensional GBM and MHD models are intrinsically nonequilibrium and 
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do not generally have an equilibrium limit: Switching off the noise cross-correlation does 
not make these models equilibrium in general ^-dimensions. In contrast, the present model 
has a well-defined equilibrium limit given by A^x = for any dimension d. Thus not only 
does our model here exhibit continuously varying universal properties, it can be driven away 
from equilibrium continuously and incrementally by tuning N^. Continuously varying uni- 
versality has been found in Ref. 19| as well. However, Ref. 19|] required coupling of the 
order parameter field with a conserved density. In contrast in our work we have the order 
parameter field only as the relevant dynamical field. Quantitative accuracy of our results 
is limited by the heuristic arguments we resorted to while evaluating the diagrams arising 
from noise cross-correlations. In order to verify this, direct numerical simulations of the 
model Langevin equations, or simulations of appropriately defined lattice-gas models should 
be performed. In the present article we have discussed the universal scaling properties at 
the critical point only. Numerical simulations of a driven 0(3) model 26| displays existence 
of spatio-temporal chaotic low-temperature regime below its critical point in the absence 
of stochastic noises. This chaos, when controlled, is replaced by spatially periodic steady 
helical states which are robust against noise. In view of these results in Ref. 26|, it would 
be interesting to examine the properties of the ordered phase below T^, and their depen- 
dences on the parameter A^x introduced above. In the above we have confined ourselves in 
discussing a usual order- disorder transition and the associated universality at the critical 
point. For any model, such a scenario holds as long as the physical dimension is greater 
than the lower critical dimension d^ of the model. For equilibrium models with continuous 
symmetries, e.g., the 0(2)-symmetric model in equilibrium, the Mermin- Wagner theorem 
tells us that di = 2. For models out of equilibrium, there are no such general results. It 
would be interesting to examine the present model, perhaps through non-perturbative meth- 
ods, at (i = 2, in particular the role and dynamics of topological defects in the presence of 
noise cross-correlations. We hope our theoretical results will inspire more detailed theoret- 
ical studies on more realistic models or experimental work on driven systems with coupled 
variables, where the role of noise cross-correlations in determining the universal properties 
may be explicitly tested. 
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